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1. Introduction 

The aim of this paper is to give some variational characterizations of the effective multiplication 
factor arising in nuclear reactor theory. This work follows a very recent paper by M. Mokhtar- 
Kharroubi [22] devoted to the leading eigenvalue of transport operators. 

In practical situations, the power distribution in a stable nuclear reactor core is determined as the 
steady-state solution of a linear transport equation for the neutron flux. Because of interactions be- 
tween neutrons and fissile isotopes, a fission chain reaction occurs in the reactor core. Precisely, when 
an atom undergoes nuclear fission, some neutrons are ejected from the reaction and subsequently shall 
interact with the surrounding medium. If more fissile fuel is present, some may be absorbed and cause 
more fissions (see for details Refs. [5, 6, 11, 25]). The linear stationary transport equation is there- 
fore of non-standard type in the sense that the source term is itself a function of the solution. When 
delayed neutrons are neglected, this equation reads 

v ■ V x cj)(x, v ) + a(x, v )(f)(x, v) — / K s (x,v,v')(j)(x,v')dfi(v') = 

Jv 

- — / K f (x,v,v')(j)(x,v')dfi(v') (1.1) 

withfree-surface boundary condition (i.e. the incoming flux is null). Here, the unknown <f>(x, v) is the 
neutron density at point i£D and velocity v € V, where V is an open subset of R^ (representing 
the reactor core) and the velocity space V is a closed subset of R , dp,{-) being a positive Radon 
measure supported by V. For the usual cases, dfi(-) is either the Lebesgue measure on R N (contin- 
uous model) or on spheres (multigroup model). The transfer cross-sections k s (-, •) and (•, •, •) 
describe respectively the pure scattering and the fission process. The nonnegative bounded function 
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cr(-, •) is the absorption cross— section [5, 11, 25]. The positive ratio k c g is called the criticality eigen- 
value (or the effective multiplication factor). It represents the average number of neutrons that go on 
to cause another fission reaction. The remaining neutrons either fail to induce fission, or never get 
absorbed and exit the system. Consequently, k c g measures the balance between the number of neu- 
trons in successive generations (where the birth event separating generations is the fission process). 
The interpretation of the effective multiplication factor k e g is related to the following three cases (see 
[5, 11,25]): 

• If k c g = 1, there is a perfect balance between production and removal of neutrons. The reactor 
is then said to be critical. 

• The reactor is sub-critical when k e g < 1. This means that the removal of neutrons (at the 
boundary or due to absorption by the surrounding media) excesses the fission process and the chain 
reaction dies out rapidly. 

• When k e ff > 1, the fission chain reaction grows without bound and the reactor is said to be 
super— critical. 

Up to now, in practical applications, the effective multiplication factor k e g was usually given by 
r a [(T — /C s ) _1 /Cy] = fceffj where the precise definition of the operators T, K s and Kf is given 
subsequently and r a [B] denotes the spectral radius of any generic bounded operator B: r a [B] = 
limj^oo Hi?™!! 1 /™. Because its requires the computation of the resolvent (T — K, s )~ l , such a char- 
acterization makes the analysis of the effective multiplication factor quite difficult to handle. In 
particular, practical estimates of k e g as well as computational approximations are rather involved and 
merely rely on (direct or inverse) power method [2]. Our main objective in this paper is to provide 
tractable variational characterizations of k c g in terms of the different data of the system and suitable 
test functions. We hope that such a characterization shall be of interest for practical computations or 
for the homogenization of the above criticality transport equation in periodic media [1, 4, 21]. We 
also believe that our characterization can be useful in the delicate optimization problem of the assem- 
bly distribution in a nuclear reactor (see the recent contribution [28] based upon the homogenization 
method and where the criticality eigenvalue k c g plays a crucial role). 

To be precise, we provide here variational characterizations of the effective multiplication factor 
of the type 



v ■ V x ip(x,v) + a(x,v)ip(x,v) — / S s (x, v , v')ip(x, v')dp,(v ') 



- — = mm ess sup 

k c S <p£W+ (x,v)GVxV 



V 



V 



v, v')ip(x, v')d/i(v') 

(1.2) 

v ■ V x ip(x, v ) + a(x, v)ip(x, v) — / H s (x,v,v')ip(x,v')dfj,(v') 



max ess inf 

, e w+ (x,v)£T>xV 



V 
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where Wp is a suitable class of positive test-functions in EPifD x V, dxdfj,(v)) (1 ^ p < oo). This 
result (Theorem 3.3) holds true under compactness assumption on the full collision operator 

JC : ij) \ — > / CE s (x,v,v') + T,f(x,v,v')) ip(x, v')dp,(v') 
Jv 

and under positivity assumptions on the fission cross-section «/(•,-,-). The main strategy to derive 
(1.2) is adapted from [22] where M. Mokhtar-Kharroubi proved similar variational characterizations 
for the leading eigenvalue of perturbed transport operators. Note that the above characterization still 
holds true for transport equations with general boundary conditions modeled by some nonnegative 
albedo operator (see Remark 3.6). 

At this point, one recalls that, besides the critical eigenvalue & e g, it is also possible to investigate 
the reactivity of the nuclear reactor core through another physical parameter, namely, the leading 
eigenvalue s(A) of the operator A = T+/C s +/C/, also associated to positive eigenf unctions. The two 
parameters k e s and s(A) are related by the following: if s(A) < then the reactor is subcritical (i.e. 
^cff < 1)> while it is super-critical whenever s(A) > 0. The reactor is critical when s(A) = 0. The 
paper [22] provides a variational characterization of the leading eigenvalue of the transport operator 
A. However, for practical calculations in nuclear engineering, the critical eigenvalue k e g is a more 
effective parameter. Actually, the existence of the leading eigenvalue s(T + K s + Kf) is not always 
ensured but is related to the size of the domain V and the possibility of small velocities (for more 
details on this disappearance phenomenon, see e.g. [19, Chapter 5]). Since the existence of the 
effective multiplication factor is not restricted by such physical constraints, it appears more efficient 
to measure the reactivity of nuclear reactor cores by k e s. This is what motivated us to generalize the 
result of [22] and provide variational characterization of k e g. 

In this paper, we also give a characterization of the criticality eigenvalue associated to the energy- 
dependent diffusion model used in nuclear reactor theory [5, 24, 25]. For this description, the critical 
problem reads 

- div x (B(x, Z)V x q(x, 0) + ct(x, Oq(x, - / E,(x, £, £') g(x, £') df 

JE 

= -j^ J V f fat,?) g(x, (1.3) 

complemented by the Dirichlet boundary conditions g(-,C)\8V = a - e - £ € E. Here E is the set of 
admissible energies £ = ^mv 2 (m being the neutron mass and v the velocity), i.e. E is a subset of 
[0, +oo[. The diffusion coefficient D(-, •) is a matrix-valued function over V x E and the unknown 
g(-, •) is nonnegative. 

The derivation of diffusion-like models for some macroscopic distribution function q(x, £) (corre- 
sponding to some angular momentum of the solution (p to (1.1)) is motivated in nuclear engineering 
by the necessity to provide simplified models tractable numerically. Such a energy-dependent diffu- 
sion model can be derived directly from a phenomenological analysis of the scattering models or it 
can be derived from the above kinetic equation (1.1) through a suitable asymptotic procedure (see 
[10] and the references therein for more details on that matter). 
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For this energy-dependent diffusion model, we give a variational characterization of fc e ff in terms 
of sup-inf and inf-sup criteria in the spirit of (1.2). 

Actually, to treat the two above problems (1.1) and (1.3) it is possible to adopt a unified mathemat- 
ical formalism. Precisely, let us denote by K s the integral operator with kernel k s (-, -, ■) and denote 
by KLf the integral operator with kernel •, •). Then, problems (1.1) and (1.3) may be written in a 

unified abstract way: 

T(f) cS + IC s (f> cS + y^-ZC/^cff = 0, (/) eS ^ 0, 

«eff 

where the unbounded operator T refers to, according to the model we adopt: 

• the transport operator: 

Tcj>(x, v) = —v ■ \7 x (fr(x, v) — a(x, v)<f>(x, v), 

associated to the absorbing boundary conditions <p\ r _ = 0. 

• the energy-dependent diffusion operator: 

Til>(x, = div x (D(x, e)V^0z, 0) ~ ZMx, 
with Dirichlet boundary conditions ip\ dvi'iO = ( a - e - £ £ E)- 

The abstract treatment of the above problem is performed in Section 2 and relies mainly on positiv- 
ity and compactness arguments. The main abstract result of this paper (Theorem 2.15) characterizes 
the criticality eigenvalue of a large class of (abstract) unbounded operators in L p -spaces. Besides this 
main analytical result, we also prove abstract results with their own interest. In particular, we pro- 
vide in Theorem 2. 12 an approximation resolution for the criticality eigenfunction 4> e g which shall 
be hopefully useful for practical numerical approximations. 

The outline of the paper is as follows. In Section 2, we describe the unified and abstract framework 
which allows us to treat in a same formalism Problems (1.1) and (1.3) with the aim of establishing 
general inf-sup and sup-inf formulae for the criticality eigenvalue of a class of unbounded operator. 
In Section 3, we are concerned with the characterization of the effective multiplication factor k e g 
associated to the transport problem (1.1). In Section 4, we investigate the effective multiplication 
factor associated to the energy-dependent diffusion model (1.3). 

2. Abstract variational characterization 

This section is devoted to several abstract variational characterizations of the criticality eigenvalue. 
It is this abstract material that shall allow us to treat in the same formalism Problems (1.1) and (1.3). 

2.1. Setting of the problem and existence result. Let us introduce the functional framework we 
shall use in the sequel. Given a measure space (ft, v) and a fixed 1 ^ p < oo, define 

x p = L p {n,dis) 

and denote by X q its dual space, i.e. X q = L q ($l,dv) (l/p + 1/q = 1). We first recall several 
definitions and facts about positive operators. Though the various concepts we shall deal with could 
be defined in general complex Banach lattices, we restrict ourselves to operators in X p (1 ^ p < oo): 
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Definition 2.1. A bounded operator B in X p is said to be irreducible if, for every nonnegative (p £ 
X p \ {0} and nonnegative ip G X q \ {0}, there exists n £ N such that 

>0, 

where (■,■) is the duality pairing between X p and the dual space X q . 
Let us denote the set of quasi-interiors elements of X p by X+ , i.e. 

X+ = {f eX p ; f(u) >0du- a.e. uo € O}. 
Notice that, if / G X+ , then (/, ip) > for any nonnegative ^ £ -Xq \ {0}. 

Definition 2.2. A bounded operator B in X p will be said to be positivity improving if its maps 
nonnegative f G X p \ {0} into Xp, i.e. 

feX p \{o},f^o^BfeX+. 

Remark 2.3. Notice that, given a bounded operator B in X p , if some power of B is positivity im- 
proving, then B is irreducible. This provides a practical criterion of irreducibility. 

Recall also several fact about power-compact operators. 

Definition 2.4. A bounded operator B in a Banach space X is said to be power-compact if there 
exists n G N such that B n is a compact operator in X. 

The following fundamental result is due to B. De Pagter [8]. 

Theorem 2.5. Let B be a bounded operator in a Banach space X. If B is irreducible and power- 
compact then r a (B) > where r a (B) denotes the spectral radius of B. 

Let T be a given densely denned unbounded operator 

T : ®{T) cXp^Xp 

such that 

s(T) < and (0 - T)-\X+) C X+. (2.1) 
Let K, s and ICf be two nonnegative bounded operators in X p . We are interested in the abstract critical 
problem: 

(T + /C s + ^-/C / )0 cfr = O, fa E 9{T) , eff > 0, fa ^ 0. (2.2) 
Since s{T) < 0, Problem (2.2) is equivalent to 

(0 - Ty l (K. s + - — )Cf)(f) e ff = fa , fa ^ 0, fa / 0. 
Let us introduce the family of operators indexed by the positive parameter 7: 

K(^)=K S + -K. f 7 >0. 

Therefore, solving (2.2) is equivalent to prove the existence (and uniqueness) of k e g > such that 1 
is an eigenvalue of (0 — T)~ 1 /C(/c e ff ) associated to a nonnegative eigenfunction. Such an existence 
and uniqueness result can be found in [19, Theorem 5.30] (see also [18]). We set /C = /C s + JCf 
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Theorem 2.6. Assume that (0 — T) _1 /C is power-compact and that (0 — T) _1 /C/ is irreducible. 
Then, the spectral problem (2.2) admits a unique solution k c g > associated with a nonnegative 
eigenf unction e g if and only if 

lim rJ(0 - T)" 1 /C(7)] > 1 and r CT [(0 - T)^ 1 ^] < 1. (2.3) 

Remark 2.7. Notice that our assumptions differs slightly from that of [19]. Actually, in [19], /? z'tf 
assumed that (0 — T) _1 /C(7) is power-compact and irreducible for any 7 > 0. Owr assumption 

implies those o/[19]. Indeed, in this case, there is an integer N G N swc/i [(0 — T) _1 /C] zj 
compact. Since, for any 7 > 0, £(7) ^ max{l, l/7}/C, owe geto fry a domination argument that 
[(0 — T) _1 /C(7)] is compact for any 7 > 0. This means that the power at which (0 — T)~ 1 ]C('y) 
becomes compact is independent 0/7 > 0. In the same way, since (0 — T)~ 1 /C(7) ^ ^(0 — T)~ l K,f 
for any 7 > 0, o«r assumption implies the irreducibility of (0 — T)~ l K(^) for any 7 > 0. Notice 
also that the result still holds if (0 — T)~ 1 /C s is irreducible. 

Remark 2.8. Under the assumptions of the previous Theorem, we point out that the mapping 7 > 
1— >• r a [(0 — T) _1 /C(7)] is continuous (see [15, Remark 3.3, p. 208]). By analyticity arguments 
(Gohberg-Shmulyan theorem), it is also strictly decreasing. Thus, fc e ff is characterized by 

r a [(0-T)- 1 lC(k e s)] = l. 

Let us now give some variational characterizations of the criticality eigenvalue k c g appearing in 
Theorem 2.6. 

2.2. Abstract variational characterization of k e g. From now on, Assumption (2.3) is assumed to 
be fulfilled. Let 

W+ :=@(T)nX+. (2.4) 

We start with the following characterization of fc e g in terms of super-solution to the spectral problem 
(2.2). 

Proposition 2.9. Assume that (0 — T)~ 1 /C(7) is power-compact and irreducible for any 7 > 0. For 

any ip G W£, let 

T + (ip) := sup{7 > such that (T + £(7)) •p is nonnegative } 
with the convention sup = 0. Then 

k ef[ = sup r+ (</?). 
<peW? 

Proof. Let 92 G X p be a nonnegative eigenfunction of (0 — T) _1 /C(fc e g) associated with the spectral 
radius r CT [(0 - T)~ l lC{k cS )} = 1, i.e. 

(o-ry 1 ic(k cS ) v = ip. 

Clearly, ip G Q)(T). Moreover, since (0 — T)~ 1 }C(k c g) is power-compact and irreducible, a well- 
known consequence of Krein-Rutman Theorem is that 

tp(u) > a.e. uj G ft 
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i.e. ip G Wp . Now 

(T + /C(fceff)) p = 

is nonnegative so that r + ((^) ^ A; e g and consequently 

A; cff ^ sup T+(ip). 

Assume now that k c g < sup r + {p) et let if) € be such that r+(^) > /c c g. Denote 7 = 

<peW+ 

T+(tp). By definition (T + fC(j))tp ^ 0, i.e. (0 - T)~ l K(-i)ip ^ if). Thus, for any n 6 N, 
[(0-T)- 1 /C( 7 )] n V ^ V'so 

r ff [(0- r)" 1 ^)] >1. 

since, fv [(0 - ^"^(fceff )] = 1, Remark 2.8 ensures that 7 ^ k c R , which contradicts the choice of 

r+(^) > k cS . □ 

The following illustrates the fact that the extremal value to the above variational result is reached 
only by the nonnegative solution to the spectral problem (2.2). 

Corollary 2.10. Assume that (0 — T)~ l K{^) is power-compact and irreducible for any 7 > 0. Then, 
for any nonnegative p € 3){T) \ {0} 

T + {p) = k e & if and only if (T + K,{k c R )) p = 0. 

Proo/ Let tp € £^(T) \ {0}, p ^ 0, be such that r+(v?) = /c cff . Then (T + /C(fc eff )) v? ^ i.e. 

(0 - T)- x 1C{k en ) p> ^ ip. (2.5) 

Suppose that 

<p + (0 - 1)- x K(k<M ) (2.6) 

and let -0* € -Xg be a nonnegative eigenfunction of the dual operator ((0 — T) _1 /C(/c e fj))* associated 
to rv[((0 - Tj^Kikes))*] = r CT [(0 - T)- 1 ^^)] = 1. Then, by (2.5) and (2.6) 

= (p, ((0 - T)- 1 /^))* r) = rv[((o - r)- 1 ^^))*] ^) 

which leads to a contradiction. Hence, 9? = (0 — T) -1 /C(/c c g) 99, i.e. 

(T + !C(k cS ))p = 0. 

Conversely, if p is a nonnegative eigenfunction of T + /C(/c c fj) associated to the null eigenvalue, then 
T + (v?) ^ fc e g and the identity T+(<p) = k e g follows from Proposition 2.9. □ 

Remark 2.11. A careful reading of the proof here above shows that, if <p £ X p is such that (0 — 
T)' l lC{k cS ) p^p^O, then p £ @(T) and (T + K(k eS )) p = 0. 
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Let us denote by e g the unique critical eigenfunction with unit norm, i.e. 4>cfi G Wp satisfies 

(T + K.(k cS ))(p c fi = (f) eS , W^esW = 1- 

Then one can prove the following approximation resolution for the criticality eigenfunction (j) e g 
whose proof is inspired by [20, Theorem 7]. Such a result shall be hopefully useful for practical 
numerical approximation of the critical mode ^> e g- of the reactor. 

Theorem 2.12. Let (ipt)k G H Xp be such that 7^ := T + (ipk) — ► k e g. We assume cp k to be 

normalized by 

\\[{Q-T)- x K(T h )] N <Pk\\ = l (keN) (2.7) 
where N is the integer given by Remark 2. 7. Let us assume that 1 is a simple eigenvalue of (0 — 
T)~ 1 /C(/c c g). Moreover, in the case p = 1, let us assume that the dual operator [(0 — T) _1 /C(/c e ff )] 
admits an eigenfunction associated to its spectral radius which is bounded away from zero. Then, 

lim \\(p k - (f> eS \\ = 

k— too 

where e g G ^ > (^") H zj f/ie unique positive eigenfunction of (0 — T) -1 /C(A; e ff ) associated to 1 
and with unit norm . 

Proof. According to the definition of 7^ := T + (ipk), (T + K,{^k))^k is nonnegative. Therefore, 

ip k ^ (0 - T)- 1 /C( 7fc )^ (2.8) 

and, iterating up to iV 

p fc < [(O-T)" 1 ^)]^^ (fcEN). (2.9) 

This shows, according to (2.7), that \\tpk\\ ^ !• Now, if 1 < p < 00 there exists a subsequence (ipk)k 
which converges weakly to some ip G X p . If p = 1, the fact that 7^ — ► fc e g combined with the 

compactness of [(0 — T)~ l K,{~fk)\ N lead to the relative compactness of 

{[(O-T)- 1 /^)^}, 

in X\. In particular, this sequence is equi-integrable and by domination (2.9), (<fik)k * s a l so equi- 
integrable. We can extract a subsequence (ipk)k converging weakly to some ip G X\. In both cases, 
the compactness of [(0 — T) _1 /C(7fc)] N together with 7^ — > fc e ff yield to the strong convergence 

[(0 - T)- 1 /c( 7fe )] * ^ -» [(0 - Ty l K,(k cS y] N ip 

so that 

|| [(O-^if^lhl, 

In particular ip ^ 0, and, taking the weak limit in (2.8), ip ^ (0 — T)~ 1 )C(k c g)ip. Now, according to 
Remark 2.11, this last inequality is actually an equality, i. e. 

Iterating again, one gets 

IHI = || [(o-r)- 1 /c(M] JV V'll = i- 
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Now, since 1 is a simple eigenvalue of (0— T)~ 1 K,(k e s), the set of eigenfunctions of (0— T)~ 1 /C(fc e ff) 
with unit norm reduces to a singleton. This shows that ijj is the (unique) weak limit of any subsequence 
of {fk)k so that the whole sequence (<pk)k converges weakly to ip G X p . The remainder of the proof 
consists in showing that the convergence actually holds in the strong sense. 

Let us consider first the case 1 < p < oo. To show now that \\ipf. — tp\\ — ► 0, it suffices to prove 
that \\(pk\\ — ► ||0||- A consequence of the weak convergence leads to 

1 = ^ liminf \\ipk\\- 

k— >oo 

Since \\ipk\\ ^ 1 f° r an Y k G N, this proves the Theorem for 1 < p < oo. 

Let us now assume p = 1. Then tpk — ► ip strongly in X\ if and only if the convergence holds in 
measure, i. e., for any Hcf! with finite ^-measure and every e > 

lim dv{oj G S ; (^(w) — tp(u)\ > e} = 0. 

Arguing by contradiction, assume there exist H C f2 with finite dz^-measure, a subsequence still 
denoted (ipk)k and some 6 > and some eo > such that 

di/{u €E;\i/> k (u)-i/>(u)\> eo}^8 for all m G N. (2.10) 

Setting 

^ fc = [(o-r)- 1 *:^)]"^ 

one has 

lim \\ip k — tp\\ = 

since ^ = [(0 - T)" 1 /C(7fc)] ^. Consequently, 

lim G 5 ; |^ fe (u;) - ip(co)\ > e /2} = 

k— >oo 

and, one deduces immediately from (2.10) that 

di/{we5;|4(w)-i/)(w)|>eo/2}^/ 2 for large A;. (2.11) 

Now, let t/;* G L°°(fi) be a positive eigenfunction of [(0 — T)~ l K,(k c R )] associated to the eigen- 
value 1, with ^* bounded away from zero. One has 



(ipk-ipk,^*) > _ (^fc(w) -V'fcM) ^*(w)di/(w) 

J{wGH; |V fe -^fc|>eo/2} 

> inf x e /2 x 5/2 = 77 > 

or else 

([(O-r)- 1 ^)]^^,^*)-^^*)^^ for large fc. 

Equivalently 

(^{[(O-r)- 1 ^)]*}^*)-^,^*)^^ for large k. 

But, {[(0 - T)- 1 A:(7 fe )]*} JV '0* converges strongly to {[(0 - T)^ 1 /C(A; cfr )]*} 7V V* in L°°(n), which, 
combined with the weak convergence tpj, — ^ ■0, implies 

{[(0 - rr^ryv) - {[(0 - Tr^Mryv) = <</>, r) ■ 
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Now, the contradiction follows from the fact that (ifi^, ip*) — > (tp, ip*) . □ 
The following characterizes k e g in terms of sub-solution to the spectral problem (2.2). 

Proposition 2.13. Assume that (0 — T)~ l K(^) is power-compact and irreducible for any 7 > 0. For 

any if 6 W£ , define 

T-(tp) :=mf{ 7 >0; -(T + %))^I p + ) 
with the convention inf = +00. Then 

k ef[ = inf r_ (</?). 

Proof. Let 99 6 be such that r_(y>) < +00, and let 7 > T_(y>). Then 

-(T + /C( 7 )) V €A-+. (2.12) 

Since (0 - T) -1 (X+) C (see Eq. (2.1)), one gets that 

<p-(0-T)- l 1C( l ) i peW+. (2.13) 

Now, let ip* G 1, be a nonnegative eigenfunction of the dual operator ((0 — T) _1 /C( 7 ))* associated 
with the spectral radius r CT [((0-T)- 1 /C(7))*] = ^[(O-T)- 1 ^)]. Note that ^[(O-T)- 1 /^)] > 
according to Theorem 2.5. From Krein-Rutman theorem, tp*(uj) > dv - a.e. uj € CI. Thus, by 
(2.13) 

(^_(0_T)- 1 /C(7)^^>>0 

i.e. 

(^V*)><(o-r)- 1 /c( 7 )^^> 

= (p, ((0 - T)- x Kb))* r) = r ff [(0 - T)- 1 /C( 7 )] . 

Since (</?, V 1 *} 7^ 0, we get r CT [(0 — T)~ 1 /C(7)] < 1. By Remark 2.8, this means that 7 > k c g. Since 
7 > T-((p) is arbitrary, we obtain 

fc e ff ^ inf T-(<p). 

¥>eW p + 

Conversely, let 7 > k cS . Then, r CT [(0 - T)~ 1 /C(7)] < 1, and 

00 

[/ — (0 — T)- 1 ^)]" 1 = £[(0 " T )" 1/C (7)] fc 

fc=0 

> [(0 - r) _1 /C(7)] n VneN. 
Given ^ e X+, set = (0 - T)" 1 ^. Then, ^ € Define 

V9= [/-(O-T)" 1 ^)]- 1 ^. (2.14) 

Clearly, <p e D(T) and 

^[(O-T)- 1 ^)]^ (neN). 
Let tp* £ X q \ {0} be nonnegative. Since (0 — T) -1 /C(7) is irreducible, there exists n S N such that 

(<p,<p*)>([(0-T)- 1 lC( r f)] n $,<p*)>0 
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Hence, (ip, ip*) > for any nonnegative p* £ X q \ {0}, so that p £ . Moreover, by (2.14), 

-(T + /C( 7 )) p = (0 — T)tp = ip£ Xp . 
Hence, 7 ^ i~-(</?) which proves that k c g ^ inf r_{p). □ 

The following result shows that only the solution of the criticality problem (2.2) realizes the above 
variational characterization. 

Corollary 2.14. Assume that (0 — T)~ l K(^) is power-compact and irreducible for any 7 > 0. Then, 
for any p E Wp 

T-((p) = k c ff if and only if (T + /C(fc e ff )) p = 0. 
Proof Let p G W+ be such that r_(p) = k cS . Then, -(T + K(k eS )) ^ so 

p^ (O-T)- 1 ^^)^. 

Arguing as in the proof of Corollary 2.10, one can prove that p = (0 — T) _1 /C(7) 99, which means 
that 

(T + K(k eB ))<p = 0. 

The sufficient condition follows directly from Proposition 2.13. □ 

We are now able to characterize the effective multiplication factor k eB by means of Inf-Sup and 
Sup-Inf criteria, where we recall that K = K s + Kf. 

Theorem 2.15. Under the assumptions of Theorem 2.6, ifKf(Xp) C Xp then the criticality eigen- 
value k e s is characterized by the following: 

1 -{T + K s )p{uj) . -(T + K t )<p(u) 
- — = mm ess sup — — - — r = max ess ml — — - — . 

K c ff <p&w+ wen <p&w+ ^en Kfp{uj) 

Proof. Let <p € W£ be given, 

T + (p) = sup{7 > ; (T + K s ) p + ^-Kf p ^ 0} = sup{7 > ; — (T + K s ) p ^ ^-Kf p} 

= sup{7 > ; — (T + K s ) p(u>) ^ — Kt p(ui) dv — a.e. u> € ft}. 

7 ' 

1 (T -\- K ) (c^) 
Since /Cf(X+) C X+, one gets — — = ess sup — r . By Proposition 2.9, 

1 . , -(T + IC s )p(cj) 
- — = mi ess sup — — r 

and the infimum is attained for the criticality eigenfunction. Similarly, let p € 

r_(p) = inf{ 7 > ; -(T + Ks) <p - -K f p > 0} = inf{ 7 > ; -(T + K s ) p > -K s p] 

= inf{7 > ; -(T + K s ) p(cu) ^ -K f p(uj) dv - a.e. uj 6 fi}. 
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So 

r_(» = inf{ 7 > ; ~ (T -± K '} f u) ^- dv - a.e. u G n}, 

1 . f -(T + /C,)v?(w) 
— — = ess ml — . 

Using Proposition 2.13, one proves that 

i i . { -(r+K>(w) 

sup — —- = sup ess ml - 



which ends the proof. □ 

2.3. The class of regular collision operators. We end this section by recalling the class of regular 
collision operators introduced in kinetic theory by M. Mokhtar-Kharroubi [19]. This class of operators 
will also be useful to study diffusion problems of type (1.3). We assume here that the measure space 
(ft, dv) writes as follows: 

ft = V x V, dv(us) = dx <8) dfx(v), u = (x, v) G ft 

where d\i is a suitable Radon measure over V . Let /C G &(L P (Q, dv)) be given by 



£:<pi — >Ktp(x,v)= / k(x,v,v')ip(x,v')dn(v') £ L p (n,dv) (2.15) 
Jv 

where the kernel k(-, •) is measurable. For almost every define 

K{x) : ^eL p (V,dfi)> — > [ k(x,v,v')ifj(v')dn(v f ) G L p (V,dfi) 



v 



and assume that the mapping /C : x G V i— > K(x) G £%(L p (V,dfi)) is strongly measurable and 
bounded, i.e. 

esssup||^(x)||^ (z ,p(v ;d/i )) < oo. 

The class of regular operators in L p spaces with 1 < p < oo is given by the following (see [19, 
Definition 4.1]). 

Definition 2.16 (Regular operator). Let 1 < p < oo. The operator K defined by (2.15) is said to be 
regular if : 

(1) For almost every x G T>, K(x) G 33(L P (V, dp,)) is a compact operator, 

(2) {fC(x) ; x G T>} is relatively compact in £$(L P (V, d[i)). 

In L 1 -spaces, the definition differs a bit. We have the following [17] 

Definition 2.17. Let K, be defined by (2.15). Then, K is said to be a regular operator whenever 
{\k(x, -,v')\ , (x, v') G T> X V} is a relatively weakly compact subset of 'L 1 (V, dp). 

The main interest of that classes of operators relies to the following (see Ref. [19] for 1 < p < oo 
and Ref. [17] for a similar result whenever p = 1): 
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Proposition 2.18 (Approximation of regular operators). Let 1 < p < oo and let K, defined by 
(2.15) be a regular operator in L P (T> x V, dx <8> dp(v)). Then, fC can be approximated in the norm 
operator by operators of the form: 

<p\ — yy~)ai(x)Pi(v) / 9i(v')(p(x,v')dfi(v') 

where I is finite, oq G L°°(V), fa G V(V, dp,) and 4 G L«(V, dp,), 1/p + 1/q = 1. 

3. The critical transport problem 

3.1. Variational characterization. This section is devoted to the determination of the effective mul- 
tiplication factor associated to the transport operator. We adopt the notations of Section 2.3, namely 
ft = V x V and du(x,v) = dx (g> dp(v). Throughout this section, we assume V to be a convex 
and bounded open subset of ~M. N while p is the Lebesgue measure over M. N or on spheres. In par- 
ticular, our results cover continuous or multi-group neutron transport problems but do not apply to 
transport problems with discrete velocities. Let 



T_ :=^(x,v)edVxV;v n(x) < 

where n(x) denotes the outward unit normal at x G d V. Let T be the unbounded absorption operator 

f T : 9{T) C X p — X p 



with domain 



ip i — ► Tip{x, v) := —v ■ V x tp{x, v) — a(x, v)ip(x, v), 
@{T) = jy> G X p ; ■ V x ip G X p and V|r_ = j. 



Here, the nonnegative function cr(-, •) G L°°(V x V) is the collision frequency. It is assumed to admit 
a positive lower bound 

a(x,v)^c>0 a.e. (x, u) G P x V. (3.1) 
Define the (full) collision operator K as the bounded linear (partial) integral operator 

/C : V ^ Xp i — > )Cip(x,v) := / Y,(x,v,v')if;(x,v')dp(v') E X p . 

Jv 

The collision kernel -, •) is assumed to be nonnegative. In nuclear reactor theory, in a fissile 
material, this collision kernel splits as 

S(x, v, v') = T, s (x, v, v') + Sf (x, v, v') 

where T, s (x,v,v') describes the pure scattering phenomena and Y,f(x,v,v') describes the fission 
processes. Define the corresponding linear operators 



K s : tp G Xp i — > )C s ^(x, v) := / S s (x, f , v')tp(x, v')dp(v') G X p 

and 

/C/ : ^ G Xp i — ► Kf4>(x,v) := / Y t f(x,v,v')ip(x,v')dfj,(v') G X p . 
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As we told it in Introduction, we are interested here in the critical problem: 

v ■ V x (p(x,v) + a(x,v)ip(x,v) — / T ls (x,v,v')(p(x,v')dfi(v') 



Jv 

= - — / Y,f(x,v,v')<p(x,v')dp(v'), (3.2) 

Kcflf Jv 

where the eigenf unction ip is nonnegative and satisfies the boundary condition <^| r „ = 0. We recall 
that the spectral bound of T is given by [29] 

s(T) = — lim inf t" 1 I a(x + sv,v)ds, 

t-*ca ( x ,v)eT>xV Jo 
t<r(x,v) 

with t(x, v) := inf{s > ; x — sv £ V}. Therefore, by (3.1), we have s(T) < 0. Moreover, 

rr(x,v) f rt ~\ 

(0 — T)~ <p(x, v ) = / exp { — &(x — vs, v)ds > (p(x — vt, v)dt, 



Jo k. Jo J 

so that (0 — T) _1 fulfills (2.1). Let us now recall the deducibility properties of (0 — T)~ l Kf for the 
continuous and multigroup models. The following result may be found in Ref. [19], Theorem 5.15, 
Theorem 5.16, (see also [29]). 

Theorem 3.1. Let T> be convex. Then, (0 — T)~ l K,f is irreducible in the two following cases: 

(1) V is a closed subset ofM. N equipped with the Lebesgue measure dp, and there exist < c\ < 
C2 < oo, such that Vq = {v 6 ; c\ < \v\ < 02} C V with 

E f (x,v,v')>0 a.e. (x,v,v') G (V x V x V ) U (V x V x V). (3.3) 

(2) V is the union ofk disjoint spheres (k ^ 1), 

k 

V = |J v u Vi = {veR N ; \v\ = n}, {n > o, i = l, . . . , k) 

1=1 

and, on each sphere, dp is the surface Lebesgue measure. Moreover, for any i, j 6 {1, . . . ,k}, 
there exists £ € {1, . . . , m} such that 

T, f (x,v,v) > a.e. (x,v,v') e (V xViXVijU (T> xVtXVj). (3.4) 

Remark 3.2. In the above case (I), corresponding to continuous models, it is possible to provide 
different criteria ensuring the irreducibility of (0 — T)~ 1 K,f (see for instance Ref. [13]). In the 
second case (2), which corresponds to multigroup transport equation, several different criteria also 
exist [23]. 

Using the notations of Section 2, we have the following characterization of the effective multipli- 
cation factor of the transport operator. 

Theorem 3.3. Let us assume that KL is a regular collision operator and that one of the hypothesis of 
Theorem 3.1 holds. The critical problem (3.2) admits a unique solution k c $ if and only if 

lim tv[(0 - r) -1 /C(7)l > 1 and r a [(0 - T) _1 /C s ] < 1. 
7^0 
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Moreover, 

v ■ V x ip(x, v ) + a(x, v)ip(x, v) — I T, s (x, v, v')ip(x, v') d/j,(v') 
- — = min ess sup 

^cff tpeWjt (x,v)£VxV 



V 



V 



v, v')ip(x, v') d/i(v') 

, (3.5) 
v ■ V x tp(x,v) + a(x,v)ip(x,v) — / T, s (x,v,v')ip(x,v') dfj,(v') 



max ess inf - 

<p£W+ (x,v)evxv 



v, v')<p(x, v') d(j,(v') 



Proof. Since K is a regular collision operator, one deduces from [19, Theorems 4.1 & 4.4] when 
1 < p < oo (respectively [17] if p = 1) that (0 — T) _1 /C is a power-compact operator in X p 
(1 ^ p < oo) under our assumptions on the measure fi. Moreover, in the continuous case, thanks to 
(3.3), for any tp £ 

Kf<p(x, v ) = / T,f(x,v,v')tp(x,v')dn(v') 
Jv 



^ / T,f(x,v,v') (p(x,v')dfi(v') > a.e. (x, v) G V x V, 
Jv 

i.e. JCf(X+) C X+. Similarly, in the multigroup case, Eq. (3.4) implies JCf(X+) C X+. Now, the 
existence of k c g follows from Theorem 2.6 while (3.5) follows from Theorem 2.15. □ 

Remark 3.4. Denote by 4> e g the nonnegative solution of (3.2), one can check that (p e ff £ W^. 
Therefore, in (3.5), the supremum and the infimum are reached for <j> = <j) e ff. 

Remark 3.5. Note that it is possible to provide practical criteria that are satisfied in nuclear reactor- 
theory and that ensure the existence ofk c g [4, 28]. Such criteria usually rely on dissipative properties 
of the pure scattering operator. 

Remark 3.6. It is important to point out that the above characterization is not restricted to the case 
of absorbing conditions but also holds for general boundary conditions modeled by some suitable 
nonnegative albedo operator. Actually, if one considers a transport operator Th associated to gen- 
eral nonnegative albedo boundary operator H which relates the incoming and outgoing fluxes in T> 
[16], then the above theorem holds true provided (0 — T^)" 1 ^ is a power-compact operator in X p 
(1 ^ p < oo) when KL S and /Cj are regular operators. This is always the case whenever 1 < p < oo 
by virtue of the velocity averaging lemma [16]. The problem is more delicate in a ^-setting and is 
related to the geometry of the domain T> [27]. 

3.2. Necessary conditions of super-criticality and sub-criticality. We shall use the result of the 
previous section to derive necessary conditions ensuring the reactor to be super-critical or sub-critical. 
Note that, for practical implications, a nuclear reactor can be operative and create energy only when 
slightly super-critical (i.e. 1 < k c g < 1 + S with S > small enough), in this case, the whole chain 
fission being controlled by rods of absorbing matter. Throughout this section, we shall assume fe e g to 
exist. 
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We shall provide lower and upper bounds on the effective multiplicative factor k e s only when the 
velocity space V is bounded away from zero. Recall that, since V is assumed to be closed, this means 
that ^ V (see also Remark 3.14). 

For almost every x € V, define K7(x) as the following operator on L P (V, dp,): 

)C(x) ^a^,^)- / f i x ' v : v ' )t ^ v '\ ^(v')dfi(v') e W{v, d/x) 

J v 1 + a{x,v)r(x,v) 

where we recall that S(x, v, v') = T* s (x, t>, v') + u, u') and t(x, v) is the stay time in V. Then, 
one defines as in [22], the following 

- . 

$ := mi ess sup -p— r 

where L p ,(V,dp) = {ip € L p (V,d[i) ; ^(«) > dp — a.e. « € V}. Then, one has the following 
estimate: 

Proposition 3.7. Under the assumptions of Theorem 3.3, ifd<\, then fc e g ^ •&. 
Proof. Assume ■& < 1. Given # G 1), let G L+(V, dp,) be such that 

ess sup T~T\ ^ ^ - 

(x,v)eVxV m\ v ) 

Let us consider then the following test-function ipo(x,v) = r(x,v)ipo(v). Since ^ V, t(-, •) is 
bounded and such an application ipo belongs to since 

t(x + tv, v) = t(x, v ) + t a.e. (x,v) E T> X V, t ^ 0, 

implies v ■ V x ipo(x, v) = tpo( v )- Then, for any 7 > 0, one sees that 

- (T + K(i))(po(x,v) (l + a(x,v)r(x,v)) (fol> (v) - [1C \x)^ Q ]{v) 



+ (1 - 0) / 7- — 7 T7 r^o (« )+ 



+ / il f w ' (3-6) 

7 ii/ 1 + <j(x,t))r(a;,f) 



Since S s ^ and 1 — # ^ 0, one sees that 

- (T + /C( 7 ))^o(x, w) ^ tf" 1 (1 + u)t(x, «)) ( ^o(^) - [K7 '(x)ifo](v) 



+ / — 7 w r^o{v)dp{v 

7 7 y l + o-(x,i;)T(a;,«) 

In particular, from the positivity of E t, one sees that, provided 7 ^ #, — (T + IC( , y))ipo(x,v) ^ for 
almost every (x, v) E T> xV. Then, from Proposition 2.13, this means that T_(ipo) ^ $ and /c c g ^ 
Since i? > 1? is arbitrary, one gets the result. □ 
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Remark 3.8. From the above result, one sees that the reactor is sub-critical whenever •& < 1. Note 
that the fact that •& < 1 implies k e s ^ 1 is already contained in [22, Theorem 7]. 

The above result provides an upper bound of k c g leading to the sub-criticality of the reactor core. 
To get a lower bound of k e g, one defines a similar quantity 

v= sup ess ml — . 

Proposition 3.9. Under the assumptions of Theorem 3.3, if$ > 1, then /c e g ^ d_. In particular, for a 
reactor core to be sub-critical, it is necessary that §_ ^ 1. 

Proof. The proof is very similar to that of Prop. 3.7. Namely, assume # > 1. For any $ G (1, ^9), let 
V>o G L p + {V,dp) be such that [/C r (£)V>o](w) ^ $il>o(v) for almost every G V x V\ Then, the 

function ipo(x,v) = T(x,v)ipo(v) belongs to and, arguing as in Prop. 3.7, thanks to Eq. (3.6), 
one sees that, (T + JC(^))ip is nonnegative for any 7 ^ Consequently, r+(v?) ^ # and Prop. 2.9 
implies that fc c g- ^ $ for any i? G (1, □ 

Remark 3.10. To the author's knowledge, the identity $ = $is an open question. Notice however 
that, according to I. Marek's result, Ref. [18], Theorem 3.2, for any x G T>, one has the identity 

. f [X7QcM(«) . [E t (sM(p) 
sup ess ml — = ml ess sup -j-^ = rJ/C (x)\, 

^L\{v,dn) v€V ^ v > ^l\{vm v&v yiy) 
where we recall that K, T {x) is an operator in LP(V, dp,). 

In the same spirit, for almost every i£D, define K T ^{x) as the following operator on L P (V, dp): 

K}( X ) :^€LP(v,d/i)^ / s/fc «y)T(sy) , , e ^ 

J 1 + cr(a;,u)T(a;,u) 

and let us define, as in [22], the set of all /3 ^ for which there exists ^ G £+(V, d/x) \ {0} such 
that 

[K7j{x)ip\(v) ^ (3i/;(v), for almost every (x, v) G V x V. 
According to [22, Lemma 4] the set / is closed so that, if one defines 

:=sup{/3,/3G/} 

then, there exists tpf G d^)\{0} such that [/Cj(x)^/](v) ^ (3fipf(v) for almost every (x, u) G 

Px7, When the velocity space is bounded away from then, j3f provides a lower bound for k e g: 

Proposition 3.11. Under the assumptions of Theorem 3.3, one has k e g /3f. 

Proof. Set <fif(x, v) = t(x, v)iftf(v) where -0/ G L P + {V, dp) is defined here above. Arguing as in the 
proof of Proposition 3.7, one sees that, since ^ V, ipf G . Therefore, Theorem 3.3 ensures that 

- — ^ ess sup j. . 

fees ( x ,v)evxv / E f (x,v,v')ip f (x,v')dp{v') 

JV 
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As in the proof of Prop. 3.7 and since Y, s and (p s are nonnegative, one gets 

1 ip f (v) 1 

— < ess sup «; — 

which ends the proof. □ 

Remark 3.12. The above Proposition provides a lower bound of the criticality eigenvalue fe e g that 
depends only on the fission collision operator ICf. In particular, a sufficient condition for the reactor 
to be super-critical is f3f > 1. 

Proposition 3.13. Let vq := infjjwj ; v & V} and let d be the diameter ofV. Define 
Af := sup essinf j Y,f(x,v,v')T(x,v')ip(v')dfi(v'). 



. uuj' 1111 If \ I 

4>eL p + (v,dti) (x,v)evxvip{v) J v 
Then, under the assumptions of Theorem 3.3, 



1 1 + (7 — 

< — H*-. (3-7) 



/c cff A f 

where a := ess sup cr(x, v). In particular, ifV bounded then 

(x,v)evxv 

k c s ^ t essinf / Hf(x,v,v')T(x,v')du(v'). (3.8) 

1 + a± (x,v)evxvj v 

Proof. Let us consider again test-functions of the form tp(x, v) = r(x, v)ip(v) where ip G d/x). 
Then, as above, according to (3.5) 

1 + a(x, v)t(x, v) - — / S s (x, u, t/)r(a;, v')ip(v')dp(v') 

1 ^ V{v)Jv 
^ ess sup zr 

k eS /--a^-.t, 



1+ ess sup a(x,v)r(x,v) 
(x,v)eVxV 



essinf / 'E,f(x,v,v')T(x,v')ip(v')dfi(v') 



(x,v)£VxV1p(v) Jy 

Since such an inequality holds for arbitrary ip(v) > and since t(x, v) ^ d/M ^ d/vo for almost 
every (x, v) ET>xV, one gets (3.7). To prove (3.8), it suffices to consider the test-function ip(v) = 1 
(v E V), which belongs to L P + {V, dp) provided V is bounded. □ 

Remark 3.14. We dealt in this section with the case of velocities bounded away from zero. For 
practical use in nuclear engineering, this is no major restriction. However, it should also be possible 
to derive explicit bounds of k c g when € V. In such a case, the exit time t(x,v) is not bounded 
anymore but behave as ^ for small \v\. Therefore, test-functions of the form ip(x,v) = r(x,v)ip(v) 

belong to if and only if^r- G £/+(V, dp,). 
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4. The critical problem FOR THE ENERGY-DEPENDENT DIFFUSION model 
4.1. Variational characterization. In this section, we are concerned with the following 

= t^- / S/(a^,0^.0^ (4-1) 

where the unknown £>(•,•) is assumed to be nonnegative and to satisfy the Dirichlet boundary condi- 
tions 

0|8d(-,O = O a.e. 

where I? is C 2 open bounded and connected subset of M. N and £7 is an interval of ]0, oof. We will 
assume throughout this section that there exist some constants o\ > (i = 1, 2) such that 

< <Ti < cr(cc,0 < CT 2 < oo, a.e. (x,£) E V x E. (4.2) 

Moreover, we assume the measurable matrix-valued application D(-, •) satisfies the following (uni- 
form) ellipticity property 

N 

essinf V dkj(x,OViVj > d i\v\ 2 (v ER N ) (4.3) 
(x,£)eVxE f-r 1 

and regularity assumption £) £ Wj^(X>) for almost every £ £ E. We will study Problem (4.1) 
in a Hilbert space setting for simplicity. Namely, set 

X 2 = L 2 (Vx E,dxd£). 

Let us assume the kernels E s (•, • , •) and E j (•,-,-) to be nonnegative and define the scattering operator 
and the fission operator 

Je 

We will assume /C s and /C/ to be bounded operators in X 2 . Define then the full collision operator 
where 

E(x, i, O = E s (x, £, O + E/(s, e, CO (x, Z,Z')€VxExE. 
Let us introduce the diffusion operator 

T : 9{T) Cl 2 ^I 2 

e 1 — ► Tq(?, = div x (B(x, (,)V x q(x, 0) ~ <t(x, £)q(x, £), 
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with domain 

®{T) = {^I 2 ; V(", g ffoP) n # 2 ( p ) a - e - C e ^ and e X 2 } 

where Hq(V) and H 2 (V) are the usual Sobolev spaces. With these notations, the spectral problem 
(4.1) reads 

(T + /C a + ^— £/)<? e ff = 0, fttf E #(T) , &fr > 0, &ff / 0. 

According to the strong maximum principle, it is clear that s(T) < and (0 — T)~ 1 (X 2 h ) C 
In order to apply Theorem 2.15, one has to make sure that (0 — T)~ l K is power-compact and that 
(0 — T)~ 1 JCf is irreducible. Let us begin with the following compactness result which is similar to 
the usual velocity averaging lemma (see [14] and [19, Chapter 2]) for transport equations and is based 
on some consequence of the Sobolev embedding Theorem [7]. 

Theorem 4.1. If K E £3{X 2 ) is regular then /C(0 — T) _1 is a compact operator in X 2 . 
Proof. By Proposition 2.18, it suffices to prove the result for a collision operator K such that 

K : geX 2 ^K Q (x,0=a(x)h(0 [ f{Z , ) e (x,Z')d£'€X 2 

J E 

where 

a E L°°(V), he L 2 (E, d£) and / E L 2 (E, df). 

Moreover, by a density argument, one can also assume / and h to be continuous functions with 
compact support in E. Let us split /C(0 — T)~ l as: 

/C(O-T)- 1 = QMiO-T)- 1 

where 

6 : q E L 2 (V,dx) i — ► [&g](x,0 = a{x)h{i)g{x) E X 2 , 
and M. is the averaging operator 

J E 

It is enough to prove that A4(0 — T)" 1 : X 2 — > L 2 (V) is compact. Let B be a bounded subset of 
X 2 . One has to show that {A4g ; g E (0 — T) _1 (i3)} is a relatively compact subset of L 2 {V). For 
any tp E B, set 

& (x,0 = (0-T)-V(a;,e). 
For almost every £ £ E, g(-, £) E Hq(T>). One extends g to the whole space M N by 



5 (x,0 ifx£V 
else. 



Clearly, for almost every £ € E, <?(■,£) E ff 1 (M 7V ). Consequently, according to [7, Proposition IX.3], 
for a. e. £ € £ and any h&R N 

\\Thg(-,0 -<?(■,£) IIl 2 (z>) < l^ll|VxK->0lli 2 (z))) 
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where T h f(x) = f(x + h) (x 6 V, h € R N ), i.e. 

( \g(x + h,0-g(x,0\ 2 dx^\h\ 2 [ \V x g(x,0\ 2 dx. (4.4) 
Jv Jv 

Now, recall that 

-&\ x (B(x,S)V x g(x,S))+a(x,€)g(x,S) =<p(x,£) (x,0 £VxE. 

Multiplying this identity by g(x, £) and integrating by paits yield, thanks to the ellipticity property 
(4.3), 



di 



\V x g{x,0\ 2 dxd^ [ \g(x,0\\<p(x,£)\dxdS. 

JVxE 



'VxE 

In particular, since B is bounded, by Cauchy-Schwarz inequality, there exists c > such that 



sup / \V x g(x,0\ 2 dxd^c. (4.5) 
-T)- 1 (B) JvxE 



ff e(o-T)-i(B) 
Then, (4.5) together with (4.4) yield 

g(x + h,0-g(x,0\ 2 dxd^c\h\ 2 . 

VxE 

By Holder's inequality, since / is continuous with compact support 

\Mg(x + h) - Mg(x)\ 2 dx = I dx 
v Jv 



<c/ \g(x + h,0-~g(x,0\ 2 dxd£ 

JVxE 



< \h\ 2 C, 

where C > does not depend on g. In particular, 

lim sup / \Mg(x + h) - Mg{x)\ 2 dx = 0. 
h ^°ge(0-T)- 1 (B) Jv 

Now, using that Mg = Mg one deduces the conclusion from Riesz-Frechet-Kolmogorov Theorem 

Ul □ 

We are now in position to prove the main result of this section where the notations of Section 2 
are adopted : 

Theorem 4.2. Let /C € SS{X<i) be regular. Assume there exists an open subset Eq C E such that 

O >0 a.e. &V x E x E . (4.6) 

Then, the problem (4.1) admits a effective multiplication factor k e g >0if, and only if, 

lim tv[(0 - r)" 1 /C(7)l > 1 and r CT [(0 - T) _1 /C s ] < 1. 
7— >o 
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Moreover, k e s is characterized 



- — = mm ess sup 

k c s -pewp (x,v)et>xe 



-div{B{x, £)V x y>(x, 0) + 0<p(x, ~ J 


E 


J 


E 



-div(B(x, £)V^(x, 0) + <K*> £M*> - J 


E 


J 


f x f (x,s;,v'Mx,e)de 

E 



(4.7) 



= max ess inf 

<p£W+ (x,t)ET>xV 

Proof. From (4.6), Kf{X^) C X% • Now, for almost every £ £ E 1 , define Tg as the following operator 
on L 2 {V): 

: ge h-> T^g(x) = div x (B{x, £)V x f>(x)) - a(x, £)e(x), 

where @(T(.) = Hq(T>) D H 2 (V) turns out to be independent of £. Since £> is connected, the 
(elliptic) maximum principle implies that (0 — Tjt) is iiTeducible (see [9, Theorem 3.3.5] or [3, 
Section 11.2]). Actually, since is the generator of a holomorphic semigroup, this implies that 
(0 — 2$) is positivity improving (see [23, p. 306]), i.e. (0 — T$) g(x) > for almost every x G T> 
provided g G L 2 (V), g(x) ^ for almost every x G T> and £ 7^ 0. Now, let ip G X2, ^>(x, £) ^ 
for almost every (x,£) e V x E, ip ^ 0. Then, /C/V ^ and (0 - T 5 )~ 1 /C / ^(x, £) > for almost 
every (x,£) £ V x E. It is easy to see that this exactly means that (0 — 1~) K.fip(x,£) > for 
almost every (x, £) G P x E and the irreducibility of (0 — T)~ l K,f follows. Since (0 — T)~ l fC is a 
compact operator by Theorem 4.1, the conclusion follows from Theorems 2.6 and 2.15. □ 

4.2. Explicit bounds. In this section, we derive explicit bounds for the effective multiplication factor 
fc e ff. As we did in Section 3.2, the strategy consists in applying Theorem 4.7 to suitable test-functions. 
We assume the hypothesis of Theorem 4.2 to be met. Moreover, we assume here that the diffusion 
coefficient D(-, •) is degenerate, i.e. 

D(x,e)=Oo(x)di(0, (x,0eVxE, 

where Oo(') is a matrix-valued application satisfying the ellipticity condition (4.3), D(-) G Wj ' (T>) 
and di(-) is a bounded real-valued application with 

ess inf di(£) > 0. 

Let Ao be the principal eigenvalue of the following elliptic problem in L 2 (D) 

fdiv x (B Or)Ve(x)) + A ^(x) = 0, (x G V) 

I ^(x) =0 (x G &D). 



(4.8) 



It is well-known [9] that Ao > and that there exists a positive eigenfunction go solution to (4.8). 
Set £+ = {tp G L 2 (E,d£,) ; > a - e - £ e In tne s P irit of Section 3.2, for almost every 
x G T>, define IC^° (x) as the following operator on L 2 (E, d£): 



7 7b Aodi(0 +v{x,Q 
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and let If be the set of all /3 ^ for which there exists ip £ £ + such that 

[k,)° ^ for almost evef y ( x ' G v x E - 

Proposition 4.3. Setting (3q := sup{/3, /3 E I}, one has k e Q ^ /3q. In particular, a necessary 
condition to the reactor to be sub-critical is (3$ < 1. 

Proof. As in Section 3.2, the set I is closed. Therefore, there exists ip G £ + such that [/C A ° (x)ip](£) ^ 
Poip(Q f° r almost every (x, £ V x E. Now, set ^o(^) C) = ^0(^)^(0; then, ip € and 

-div x (B(x,0^Mx,0) = -di(£)il>(£)aiy x (Po(x)VQo(x)) 

= Aodi(£M£)eo(aO (x,Z)eVxE. 

Consequently, thanks to (4.7) one has 



1 

- — ess sup 

< ess sup — 7 . = ess sup 



[Ao<M£) +(t(x,£)]qo(x)iJj(£) - Qq(x) j 


f X s (x,£,OV>(£'R' 

E 


Qo(x) j 


E 



ii 



which proves that ^ ^ ^ . □ 

In the same spirit, for almost every x € V, define /C A °(x) as the operator on L 2 (E, d£) given by 

/C A °(x) : ^ G L 2 (E, d£) i— > / - WW G L 2 (E,<%) 

Je Xod^) + a{x,£) 

where we set S(cc, £, £') = £, £') + £, £'). As in Section 3.2, set 

inf esssup J^»l, md , := sup e88 inf I^MffiQ. 

Then, one has the following bounds of k c g, in the spirit of Propositions 3.7 & 3.9. 

Proposition 4.4. Under the assumptions of Theorem 4.2, if$> 1, then k e $ ^ 0« the other hand, 
if$> 1, then k cS ^ 

Proof. The proof is very similar to that of Prop. 3.7 & 3.9. We only prove the first part of the result, 
the second part proceeding along the same lines. Assume thus that ?? > 1. For any i9 € (1,^), let 

ipo £ £ + be such that ess inf ^ ^)^x>xE ~ °io(g)°^ ^ ^' Choose then the test-function £) = 
^0(^)^0(0- Such an application </? belongs to and, as in the above proof, for any 7 > 0, one 
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sees that 

(T + JC( 7 )M*, = ^ (Aodi(e) + (x, 0) ( - ^o(0 + [/C Ao 

(0 - 1)[/C S A °(*)M£) + ^[/C A °(x)MO 

where /C A °(:r) has been already defined and the definition of Kg°(x) is similar (S s replacing E^). 
Then, from the positivity of S s and S the assumption $ > 1 implies that (T+/C(7))<^ is nonnegative 
for any 7 ^ 1?. Consequently, r+(<£>) ^ $ and Prop. 2.9 implies that k c g ^ Since i9 € (1,#) is 
arbitrary, one obtains fc e g ^ □ 

Whenever E is of finite Lebesgue measure, one has the following practical criteria, already stated 
by C. V. Pao [24, Theorem 5.3] using completely different arguments. 

Corollary 4.5. Assume E to be of finite Lebesgue measure. If 

Ao di(£) + a(x, < [ [£ s (z, £, O + £, £')]< O^xE, (4.9) 
f/ierc, the reactor core is non super-critical, i.e. k e $ ^ 1. 

Proof. Since E is of finite Lebesgue measure, the constant function ip = 1 E such that = 1 f° r 

any ( £ £ belongs to £ + . Then, assumption (4.9) means exactly that [K. x ° (x)l E ](0 > for 

almost any (x, £ V x E. Therefore, $ > 1 and the conclusion follows from Prop. 4.4. □ 

Remark 4.6. Notice that, under the above assumption, one has 

, . . , f E [z s (x,tj,e) + x f (x,q,ow 

fceff ^ ess ml — — — ; — — . 

(*,{) Aodi(0 + ^.0 
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